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Abstract —Silencer was considered with and without 
modifying in order to reduce the vibration. Design and 
modelinghave been done with specifying different material 
properties. Finite Element Method was used for the 
modeling and simulations. The harmonic analysis has been 
performed by using ANSYS 18.0. The natural frequency 
and working frequency are the very important parameters 
to study the resonance. It is mandatory to avoid this 
resonating condition. These frequencies are distinguished 
with the help of ANSYS 18.0. 
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I. INTRODUCTION 

Silencer looks like a cantilever structures which is one of 
the important parts of the exhaust system. At running 
condition of vehicle Silencers and Mufflers are subjected 
to various structural, thermal and vibration loads. The 
Silencer may get afailure or damaged because of many 
reasons, vibration from engine and road excitations are 
amongst them. The vibration failure happens mainly 
because of the resonant frequencies occurring in the 
defined frequency range. If we consider the road 
excitation, vertical accelerations are dominant in Mufflers 
[7]. The selection ofexhaust silencer is depends on the 
level of reduction in the engine exhaust noise. In order to 
select the exhaust silencer, we have to compromise 
between the mechanical, aerodynamic, structural and 
acoustical performance in conjunction with the cost of the 
resulting system. Hence, it is more important study the 
vibrations in deep which would further help to improve the 
efficiency of Silencer, improving life, minimize cracks and 
life of a silencer [lj.The main objective to achieve 
Vibration optimization of Silencer using FEA 
(ANSYS).Silencers and Mufflers are used to control the 
engine exhaust noise. Mufflers and Silencers work in a 
similar way. The technical functionality between these two 
is very similar. A silencer is a traditional name given to a 
device which attenuates noise, Muffler is a device, 
designed to minimize the engine exhaust noise [3]. The 
basic function of the silencer is to reduce noise and pass 
exhaust generated in the engine 


II. NEED FOR ANALYSIS 

In India, the Automobile silencer under study belongs to a 
popular 2-Wheeler manufacturer with the rated HP of the 
engine up to @7.69HP. at very high speed and 
temperature exhaust gases coming out from the engine. 
Silencer has to reduce noise, vibrations. While doing so it 
is subjected to thermal, vibration and fatigue failures 
which cause cracks. So it is necessary to analyze the 
vibrations which would further help to pursue future 
projects to minimize cracks, improving life and efficiency 
of the silencer. [1] 

III. LITERATURE REVIEW 

The first stage of design and analysis of exhaust system 
had done [1]. For the modeling of the exhaust system, he 
used conventional FEM package. The material property 
was the main parameter to study. To avoid the resonating 
condition, FFT analysis was carried out to distinguished 
natural frequency and working frequency. From his study 
it is cleared that to increase the dynamic performance, the 
thickness of the different parts plays an important role [1]. 
The advantages of different designs of Mufflers had listed 
[7]. He took into consideration different functional 
requirements of muffler like style, shape, size, cost, 
desired sound, durability, backpressure and adequate 
insertion loss [7]. 

The experimental and numerical study had performed for 
the tailpipe noise of a muffler for a wide range of throttle 
acceleration [6]. Thevalidation was done for the transient 
acoustic characteristics of a muffler in theanechoic 
chamber according to the Japanese Standard (JIS D 1616). 
At the 2 nd order of engine rotational frequency, the 
simulation results are in good agreement with the 
experimental data CFD simulations had performed to 
calculate transfer matrix of an engine exhaust muffler [5]. 
The study considers the effect of with an without mean 
flow on the acoustic performance of a silencer. The 
developed CFD model has theability to predict the acoustic 
performance of a complex muffler. It also has theeffecton 
heat transfer to produce reasonable results of the exhaust 
noise [5]. 

Next study considers weak shock wave propagation inside 
silencer of exhaust pipe [4]. To solve Euler, unsteady, 


www.ijaers.com 


Page | 5 





International Journal of Advanced Engineering Research and Science (IJAERS) [Vol-5, Issue-1, Jan- 2018] 

https://dx.doi.Org/10.22161/iiaers.5.l.2 ISSN: 2349-6495(P) / 2456-1908(0) 


compressible, two-dimensional equations, 2nd order total 
variation diminishing (TVD) scheme was used. To study 
the effect of silencer configuration on weak shock wave 
propagation eight different models of the silencer used. At 
the inlet incident plane shock was assumed and its Mach 
no. was varied from 1.01 to 1.30. The pressure distribution 
and velocity field of flow system wereanalyzed on 
designed silencer [4]. 

IV. MODAL ANALYSIS OF EXISTING MODEL 

Silencer model has been taken for the modal analysis. 
Table 1 shows the some of the properties which have been 
exclusively used in the simulation. Figure 1 shows the 
geometry file which has been drawn and model by using 
Ansys package. Figure 2 shows the meshing of existing 
silencer model. For the meshing Tetrahedral elements were 
used. The existing silencer model is somewhat complex, 
and because of this complexity, a simple tetrahedral mesh 
element was selected while meshing in Ansys. In the 
numerical simulations, geometry and meshing are the very 
crucial steps. So, special attention was given towards the 
modeling of thesilencer. After these two main steps, 
understanding of the physics is mandatory. For solving this 
complex system, initial and boundary condition plays an 
important role. Figure 3. shows the boundary condition 
used in the simulations. After meshing,simulations have 
been performed to obtain the model analysis for first five 
natural frequencies. These frequencies are shown in Table 
2. The contours of natural frequencies of theexisting model 
have been shown in Figure. (4-8). On comparing all these 
natural frequencies, maximum bending was observed in 
mode 4. 


Table 1.Material Properties. 


Element type 

Tetrahedral 

Material 

Steel 

Mesh type 

Solid Mesh 

Density 

7.85e-6 Kg/Mm3 

Poisons ratio 

0.3 

Yield Strength 

520Mpa 




Fig. 2.Meshing of silencer mode 



Q: ORIGINAL JEPIPE 

Fixed Support 
Time: 1. s 


| | Fixed Support 


F =0 


Fig. 3. Boundary conditions for existing silencer model 


Geometry and meshing were done using ANSYS 18.0. 
After meshing simulations has been performed to obtain 
the model analysis for first five natural frequencies. These 
frequencies are shown inTable 2. 


Table 2.Natural Frequency at first 5 modes for anExisting 
silencer. 


Mode 

1st 

2nd 

3rd 

4th 

5th 

Frequency 

(Hz) 

146.58 

287.71 

294.96 

435.74 

574.54 


R: ORIGINAL_Modal 



Total Deformation 
Type: Total Deformation 
Frequency: 146.58 Hz 
Unit: mm 


20.604 Max 

18.314 

16.025 

13.736 

11.447 

9.1572 

6.8679 

4.5786 

2.2893 


OMin 


Fig. 4. mode frequency for 146.58Hz 


Fig. 1. Silencer existing model. 
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R: ORIGINAL Modal 


Total Deformation 2 
Type: Total Deformation 
Frequency: 287.71 Hz 
Unit: mm 

■ 53.486 Max 

47.543 
— 41.6 
— 35.657 
— 29.714 
-• 23.771 
— 17.829 
— 11.886 
m 5.9428 



0 Min 


Fig. 5. mode frequency for 287.71Hz. 


R: ORIGINALModal 

Total Deformation 3 
Type: Total Deformation 
Frequency: 294.96 Hz 


■ 48.467 

— 42.408 

— 36.35 

— 30.292 

— 24.233 

— 18.175 

— 12.117 

— 6.0583 
B OMin 



Fig. 6.3 rd mode frequency for 294.96Hz. 


R: ORIGINAL Modal 



Total Deformation 4 
Type: Total Deformation 
Frequency: 435.74 Hz 
Unit: mm 


17.832 Max 

15.851 
13.87 
11.888 
—I 9.9068 
— 7.9255 
—I 5.9441 
3.9627 
1.9814 


0 Min 


Fig. 7.4 th mode frequency for 435.74Hz 


R: ORIGINAL_Modal 



Total Deformation 5 
Type: Total Deformation 
Frequency: 574.54 Hz 
Unit: mm 


17.652 Max 

15.69 
13.729 
—I 11.768 

— 9.8064 

— 7.8451 
—| 5.8838 

3.9226 
1.9613 


OMin 


Fig. 8.5 th mode frequency for 574.54Hz. 


has been applied to the model. The boundary condition 
used in the simulation is shown in figure 9. In 
Table 3, boundary conditions are specified. 



S: ORIGINAL JHarmonic Response 

Harmonic Response 
Frequency: 0. Hz 


| Acceleration: (Real) -49.05, (Imag) 0. m/s z 
Components: (Real) -49.05 ,0. ,0. m/s 2 
Components: (Imag) 0. ,0. ,0. m/s 2 


Fig. 9. Boundary condition of existing model 


Table 3. Boundary condition of existing model. 


Frequency spacing 

Linear 

Range Minimum 

0 Hz 

Range Maximum 

1000 Hz 

Solution Interval 

10 

User Defined Frequencies 

Off 

Solution method 

Mode Superposition 

Include Residual Vector 

No 

Cluster Result 

No 


Figure 10. gives the clear idea regarding the amplitude in 
Mega-Pascal (MPa) and natural frequencies of different 
modes. The maximum frequency was observed in mode 
four, Hence it corresponds to the maximum amplitude of 
13.981 MPa. For the same mode, the equivalent stress is 
found 2.143 MPa (See Figure 11). 


S: ORIGINAL_Harmonic Response 

Equivalent Stress 

Type: Equivalent (von-Mises) Stress 
Frequency: 1000. Hz 
Sweeping Phase: 0. “ 

Unit: MPa 


2.143 Max 

1.905 
1.6669 
1.4288 
1.1907 
0.95265 
0.71457 
0.47649 
0.23841 

0.00033004 Min 



Fig. 10.Equivalent Stress of existing model 


V. HARMONIC ANALYSIS OF EXISTING MODEL 

The harmonic analysis has been performed on the existing 
model. The mode 4 was selected for the harmonic analysis 
because it shows maximum bending. The excitation of 5G 
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Frequency (Hz) 

Fig. 11. Acceleration frequency response of the existing 
model. 



Fig. 12.Maximum principal stress response of existing 
model 

VI. NEED OF MODIFICATION OF SILENCER 

According to Jis D 1601 (Japanesestandard), 33 Hz to 67 
Hz is the damageable frequencies of vibration testing for 
automobile silencer. Vibration frequencies are necessary to 
be reduced, these causes more vibration hence noise in the 
exhaust system. Mode 4 shows more bending, so we took 
mode 4,an optimized modal for reducing the vibration. For 
reducing the failure in the existing model we add single 
stiffener in the existing model. 

VII. ADDITION OF SINGLE STIFFENER IN 
OPTIMIZED MODEL 

Single Stiffener is added to shift the frequency of the 
optimized model. First, we add a single stiffener and 
analysis results. The analysis was done and first five 
natural frequencies are determined and tabulated in Table 
4. Figurel3. Shows the modified geometry in which a 
stiffener is added. After modifying the geometry, thenew 
mesh was generated with the Ansys package. Fig.l4shows 
the meshing of the modified geometry. The boundary 
condition for the modified model is shown in Fig. 15. 



Fig. 13.Modified Model by adding a single stiffener plate 



Fig. 14. Discretized Modified Model by adding a single 
stiffener plate 


N: SINGLE STEFFNERJEPIPE 

Fixed Support 
Time: 1. s 

Fixed Support 



Fig. 15. Boundary conditions for modified model by 
adding a single stiffener plate 

VIII. MODAL ANALYSIS OF OPTIMIZED MODEL 

The modification was done in the existing silencer model 
with adding a stiffener plate. After this modification, 
modal analysis has been performed and itsresults in ten 
natural frequencies. The contour of total deformation is 
shown for different natural frequencies in Figure. (16-25). 
The first ten natural frequencies after the modification is 
tabulated in Table 4. 


www.ijaers.com 


Page | 8 












International Journal of Advanced Engineering Research and Science (IJAERS) 
httos://dx. doi. org/10.22161/ijgers. 5.1.2 


[Vol-5, Issue-1, Jan- 2018] 
ISSN: 2349-6495(P) / 2456-1908(0) 


O: SINGLE STIFFNER .Modal 



Total Deformation 
Type: Total Deformation 
Frequency: 146.81Hz 
Unit: mm 


20.59 Max 

18.302 

16.014 

13.727 

11.439 

9.1511 

6.8634 

4.5756 

2.2878 


0 Min 



O: SINGLE STIFFNER _Modal 

Total Deformation 5 
Type: Total Deformation 
Frequency: 587.02 Hz 
Unit: mm 


18.01 Max 

16.009 
14.008 
12.007 
10.006 
8.0046 
6.0034 
4.0023 
2.0011 


OMin 


Fig. 16.1 st mode of modified silencer 


Fig. 20.5 th mode of themodified silencer. 


O: SINGLE STIFFNER _Modal 

Total Deformation 2 



Type: Total Deformation 
Frequency: 287.74 Hz 
Unit: mm 


53.489 Max 

47.546 

41.603 

35.66 

29.716 

23.773 

17.83 

11.887 

5.9433 


OMin 


Fig.l7.2 nd mode of themodified silencer. 


O: SINGLE STIFFNER .Modal 

Total Deformation 3 



Fig. 18.3 rd mode of themodified silencer. 


O: SINGLE STIFFNER .Modal 

Total Deformation 4 



Type: Total Deformation 
Frequency: 439.13 Hz 
Unit: mm 


17.685 Max 

15.72 
13.755 
11.79 
9.8252 
7.8602 
5.8951 
3.9301 
1.965 


OMin 


Fig. 19.4 th mode of themodified silencer. 


O: SINGLE STIFFNER .Modal 



Fig.21.6 th mode of modified silencer by adding single 


stiffener. 


O: SINGLE STIFFNER .Modal 



Total Deformation 7 
Type: Total Deformation 
Frequency: 781.75 Hz 
Unit: mm 


50.868 Max 

45.216 
39.564 
33.912 
28.26 
22.608 
16.956 
11.304 
5.652 


0 Min 


Fig. 22. 7 th mode of modified silencer by adding asingle 


stiffener 



O: SINGLE STIFFNER .Modal 

Total Deformation 8 
Type: Total Deformation 
Frequency: 953.09 Hz 
Unit: mm 


3.8003 

0 Min 


Fig.23. 8 th mode of modified silencer by adding single 
stiffener 


34.202 Max 

30.402 
26.602 
22.802 
19.001 
15.201 
11.401 
7.6005 
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O: SINGLE STIFFNER .Modal 


Total Deformation 9 




Max 


Type: Total Deformation 
Frequency: L067. Hz 
Unit: mm 

■ 40.083 Max 

35.629 


31.176 
26.722 
22.268 
17.815 
13.361 
8.9073 
4.4537 
0 Min 


Fig. 24. 9 th mode of modified silencer by adding asingle 


stiffener 


O: SINGLE STIFFNER .Modal 



Total Deformation 10 
Type: Total Deformation 
Frequency: 1403. Hz 
Unit: mm 


52.452 Max 

46.624 
40.796 
34.968 
29.14 
23.312 
17.484 
11.656 
5.828 


0 Min 


Fig.25. 10 th mode of modified silencer by adding single 


stiffener. 


Table 4.Natural Frequency at first 10 modes for 
anoptimized silencer. 


Mode 

1 st 

2 n 

d 

3 r 

d 

4 

t 

5 

t 

6 

t 

7 

t 

8 

t 

9th 

10 

th 

Freque 

14 

28 

29 

4 

5 

7 

7 

9 

10 

14 

ncy 

6. 

7. 

4. 

3 

8 

7 

8 

5 

67 

03 


IX. HARMONIC ANALYSIS OF OPTIMIZED 
MODEL 

After the adding a single stiffener plate again the harmonic 
analysis has been done in thesame manner as previously 
done for the existing model. Table 5 shows the boundary 
condition used in the modified model. The same procedure 
has been performed for the harmonic analysis. The 
excitation of 5g had applied to the modified model. Figure 
26 shows the boundary condition used in the modified 
model. 



P: SINGLE STIFFNER .Harmonic Response 

Harmonic Response 
Frequency: 0. Hz 


| Acceleration: (Real) -49.05, (Imag) 0. m/s z 
Components: (Real) -49.05 ,0. ,0. m/s z 
Components: (Imag) 0. ,0. ,0. m/s 


Fig. 26.Boundary Conditions of optimized model. 


Amplitude vs. frequency plot for the modified model after 
the addition of single stiffener plate is shown in Figure 29. 
The maximum amplitude of the stress is found 7.0695 
MPa at 435 frequency. The corresponding equivalent 
stress is 1.9199 MPa (see Figure 29). The acceleration and 
frequency response is also plotted in Figure 27. The 
maximum acceleration of 156.85 m/s 2 is obtained at the 
natural frequency of 435 Hz. 



Fig.27.Acceleration frequency response of the optimized 


model. 



P: SINGLE STIFFNER .Harmonic Response 

Equivalent Stress 

Type: Equivalent (von-Mises) Stress 
Frequency: 1000. Hz 
Sweeping Phase: 0. ° 

Unit: MPa 


1.9199 Max 

1.7066 
1.4933 
1.28 
1.0667 
0.85345 
0.64016 
0.42688 
0.21359 
0.00030285 Min 


Fig. 28.Maximum principal stress response of optimized 


model 
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Frequency (Hz) 


Fig. 29.Maximum principal stress response of optimized 
model 

X. EXPERIMENTAL VALIDATION 

The silencer model has been taken for the experiment. The 
experiment includes the testing of vibration by 
theaccelerometer.Fig.30 shows the experimental setup for 
amodified model having single stiffener. The sensor was 
mounted at the top of the modified silencer as shown in 
Fig.30. To impose the vibration to the model a hammer has 
been applied to the silencer. Because of the hammering 
action, the vibrations get induced in the model. These 
vibrations are then sensed with the sensor and 
corresponding readings are shown by theaccelerometer. 


Table 5.First six modal frequency of vibration by FFT 
analyzer 


Mode 

1 

2 

3 

4 

5 

6 

(Frequency 
in Hz) 

97 

391 

614 

775 

958 

1066 


Frequency response function of the exhaust system is 
shown in Fig.31for different modes. 



Ho OlW 0iH> CK3.H5 

feroisar llfiO.156 narfL007^s A 2 racOJJH m/s A 2 msr:JJ05 ry's^ 

Fig. 31. Frequency response function of the exhaust system 
for different modes. 



Fig. 30.Experimentation System 


The experimental validation is done by using FFT (Fast 
Fourier Transform) analyzer. The FFT spectrum analyzer 
samples the input signal computes the magnitude of its 
sine and cosine components and displays the spectrum of 
these measured frequency components. The advantage of 
this technique is its speed. Because FFT spectrum 
analyzers measure all frequency components at the same 
time, the technique offers the possibility of being hundreds 
of times faster than traditional analog spectrum analyzers. 
The result obtained by FFT analyzer for first six natural 
frequencies are determined and tabulated as follow. 


XI. RESULT 

Table 17 shows the comparison of the natural frequencies 
of vibration silencer by FEM package and FFT analyzer. 
The comparison shows that the natural frequency by both 
methods is nearly same (see Fig.32). 


Table 6. Six modal frequency of vibration. 


Mode 

FEM (Frequency 
in Hz) 

FFT(Frequency in 
Hz) 

1 

146 

97 

2 

439 

391 

3 

587 

614 

4 

770 

775 

5 

953 

958 

6 

1067 

1066 



Mode 

Fig.32. Comparison ofFEA and testing 
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XII. CONCLUSION 

The silencer natural frequencies have been calculated by 
using the ANSYS package and by FFT analyzer. By both 
the method the natural frequencies are nearly same and 
that are useful while the design of silencer to avoid the 
resonance. The natural frequency silencer has been 
calculated by using the Ansys 18.0. The modification in 
the design has enhanced the dynamic performance of a 
silencer. Particular modification in design incorporates the 
addition of stiffener in the modified silencer. The 
amplitude of acceleration is reduced from 289.28 m/s 2 to 
156.85 m/s 2 due to theaddition of stiffeners. Maximum 
Principal Stress value reducedto alower value than existing 
design.from 13.981Mpa to 7.0695Mpa. Vibration modes 
of silencer haven't been drastically altered, but stiffener 
method can be used to reduce excitations and 
helptoincrease fatigue life and efficiency of the silencer. 
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